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Microorganisms often encounter anisotropy, for example in mucus and biofilms. We study how anisotropy and elasticity of 
the ambient fiuid affects the speed of a swimming microorganism with a prescribed stroke. Motivated by recent experiments 
on swimming bacteria in anisotropic environments, we extend a classical model for swimming microorganisms, the Taylor 
swimming sheet, actuated by small-amplitude traveling waves in a three-dimensional nematic liquid crystal without twist. We 
calculate the swimming speed and entrained volumetric flux as a function of the swimmer’s stroke properties as well as the elastic 
and rheological properties of the liquid crystal. These results are then compared to previous results on an analogous swimmer 
in a hexatic liquid crystal, indicating large differences in the cases of small Ericksen number and in a nematic fiuid when the 
tumbling parameter is near the transition to a shear-aligning nematic. We also propose a novel method of swimming in a nematic 
fiuid by passing a traveling wave of director oscillation along a rigid wall. 


1 Introduction 

The nature of the fiuid through which a microorganism swims 
has a profound effect on strategies for locomotion. At the 
small scale of a bacterial cell, inertia is unimportant and 
locomotion is constrained by the physics of low-Reynolds- 
numbei[^fiowsM3. In a Newtonian liquid such as water, low- 
Reynolds number locomotion is characterized by two distinc¬ 
tive properties: a vanishingly small timescale for the diffusion 
of velocity, and drag anisotropy, which is a difference between 
the viscous drag per unit length on a thin filament translating 
along its long axis and transverse to its long axis® In resistive 
force theory, drag anisotropy is required for locomotion®®. 

In complex fluids such as polymer solutions and gels, the 
elasticity of the polymers introduces a new timescale, the 
elastic relaxation timescale, which is much longer than the 
timescale for the diffusion of velocity®. When the fluid has 
an elastic response to deformation, swimming speeds can in¬ 
crease or decrease dependingon the body geometry and the 
elastic relaxation ti mescal e®tBI^ and the so-called scallop the¬ 
orem does not applyC^Gll. Swimmers can move faster in gels 
and networks of obstacles than in a Newtonian liquid 
When the flagellum size is similar to the size of the polymers, 
local shear-thinning may be the mimary cause of swimming 
speed variations in such fluidsE2JE3 
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1 The Reynolds number Re for the flow of a fluid with viscosity p, density p, 
characteristic flow length L, and characteristic flow velocity v is Re = pvL/p. 


Like polymer solutions and gels, liquid crystals have an 
elastic relaxation time scale, but they also alter the drag 
anisotropy required for propulsion since the fluid itself ex¬ 
hibits anisotropy. For example the nematic liquid crystal phase 
consists of rod-like molecules which spontaneously align in 
the absence of an external field. The consequences of molec¬ 
ular anisotropy on the locomotion of microorganisms have re¬ 
cently been explored experimentally. Proteus mirabilis cells 
were found to align wit h the nem atic director field and form 
multi-cellular assembliesE^l^^ini (pig^ la). When swimming 
near nematic droplets, surface topological defects were shown 
to play an important role in bacterial escape from the liquid 
crystal interfacel^. Collective dynamic effects and director- 
guided motion was also observed in Bacillus suhtilis at low 
bacterial volume fraction, and a local melting of the liquid 
crystal caused by the bacteria was found^ (Fig. lb,c). Po¬ 
tential applications include the delivery of small cargo us¬ 
ing the direction of molecular orientational. Understanding 
these results may be relevant in understanding locomotion in 
biofilmsl^(Fig. Id), and is complementary to recent work on 
active nematics, or soft active matter, in which dense suspen¬ 
sions of microorganisms themselves can exhibit nematic-like 

orderin^SHSl (Pig le). 

A classical mathematical model of swimming microorgan¬ 
isms is Taylor’s swimming sheet®, in which either transverse 
or longitudinal waves of small amplitude propagate along an 
immersed sheet of infinite extent. Extensions of this model 
have been used to study other irnpor tant phenomena such as 
hydrodynamic syn chroni zation^^^^M] interactions with other 
immersed structures^^^l^ and geometric optimizationSH. Other 
variations on this asymptotic model have been used to study 
locomotion in a wide variety of complex fluids by numerous 
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Fig. 1 (Color online) (a) Dynamic assembly of bacteria in disodium cromoglycateE^ (b) Melting of the liquid crystal medium behind a 
flagellated swimmeiG3. (c) Disclinated texture observed as a collection of bacteria locomote in a nematic liquid crystall^. (d) Aligned 
Pseudomonas aeruginosa cells in the liquid crystal matrix of concentrated dnaES. (e) Aligned swimmers in a bacterial flock. Figure from 
Gregory Velicer (Indiana University Bloomington) and Juergen Bergen (Max-Planck Institute for Developmental Biology), (f) Sketch of a 
swimming sheet (not to scale) immersed in a nematic liquid crystal with director held n(v,y,t). The propagating wave has wavelength iTzjq, 
small amplitude b <^2%/q, and wave speed c = (o/q. The director held n makes an angle 0 with the x axis. 


authors®. Locomotion in liquid crystals, however, h as not 
yet seen much theoretical treatment. In previous works®®, 
we studied a one-dimensional version of Taylor’s swimming 
sheet in a two-dimensional hexatic LC film. Departure from 
isotropic behavior in that model is greatest for large rotational 
viscosity and strong anchoring boundary conditions, and the 
swimming direction depends on fluid properties. Further un¬ 
usual properties for Taylor’s swimming sheet were observed, 
such as the presence of a net volumetric flux. Because the 
nematic phase is more commonly observed than the hexatic, 
the present study is intended to explore new features that arise 
with nematic order, and also to determine when, if ever, the 
hexatic model can be used to accurately describe swimming 
in a nematic liquid crystal. 

In this article we extend the Taylor swimming sheet model 
to the study of force- and torque-free undulatory locomotion 
in a three-dimensional nematic liquid crystal, with tangential 
anchoring of arbitrary strength on the surface of the swim¬ 
mer. We assume the director lies in the xy-plane and does not 
twist (Fig. [T). Alternatively the problem could be considered 
as filament motion in a two-dimensional nematic fluid. By 
performing an asymptotic calculation to second-order in the 
wave amplitude, assumed small compared to the wavelength. 


we examine how fluid anisotropy and relaxation affects swim¬ 
ming speed. We show how the swimming velocity depends on 
numerous physical parameters, such as the rotational viscos¬ 
ity 7, anisotropic viscosities /i/, the Frank elastic constants Ki, 
the tumbling parameter A, and the Ericksen number Er, which 
measures the relative viscous and elastic forces in the fluid. 
The rate of fiuid transport induced by swimming is also inves¬ 
tigated; unlike in a Newtonian fiuid, the induced fiuid fiux can 
be either along or against the motion of the swimmer. 


The paper is organized as follows: In §2.1| we describe the 
stresses that arise in a continuum treatment of a nematic liquid 
crystal near equilibrium. In §2.2| we use these stresses to de¬ 
rive a set of coupled equations for the flow field and local ne¬ 
matic orientation. Eollowing Tayloi^, we nondimensionalize 
and expand these equations perturbatively to first- and second- 
order in wave amplitude and derive an integral relation for the 
swimming speed and volume fiux in §2.3| and §2.4[ The de¬ 


pendence of the swimming speed and fiux on Ericksen num¬ 
ber, rotational viscosity, and tumbling parameter is described 
in ^ In ^ ^ we show that a propagating wave of director os¬ 
cillation can result in fiuid pumping and locomotion of a pas¬ 
sive fiat surface. To determine the regimes in which the results 
for swimming speed and fiux are comparable in nematic and 
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hexatic fluids, and where they differ, we plot these quantities 
side-by-side and discuss the results in the Discussion, Q 

2 Theory 

2.1 Viscous and elastic stresses 

In a continuum treatment of a nematic liquid crystal, a local 
average of molecular orientations is described by the director 
field n. The fluid’s viscous stress response to deformation is 
approximated by incorporating terms linear in the strain rate 
that preserve n ^ — n symmetry. In a n incompressible ne¬ 
matic, the deviatoric viscous stress®® is 

= 2/iE + 2/iinn(n- E -n) (nE -n + n- En), (1) 

with E = [Vv + (Vv)^] /2 the symmetric rate-of-strain tensor, 
and V the velocity field. The shear viscosity of an isotropic 
phase is jl, and and are viscosities arising from the 
anisotropy. The coefficients and can be negative, but 
the physical requirement that the power dissipation be positive 
yields bounds of /i > 0, Ail > -IpL, and Ail* + Ml > -3Ai/2. 
While both the nematic phase and the hexatic phase we studied 
previously® are anisotropic, the hexatic phase has = 

0 and thus has an isotropic viscous stress tensor, in contrast 
with the nematic. 

Meanwhile, the elastic free energy for a nematic liquid crys¬ 
tal is 

^=^(V.n)2+^(n-Vxn)2+^[nx(Vxn)f, (2) 

where K\ is the splay elastic constant, K 2 is the tw ist elastic 
constant, and ^3 is the bend elastic constant®®. The total 
free energy in the fluid (per unit length) is Fqi = f ^dxdy. 
As mentioned earlier, for simplicity we do not consider twist, 
and thus we disregard K 2 . Thus, the angle field 6{x,y^t) com¬ 
pletely determines the nematic configuration (Fig. If). Com¬ 
paring again with our previous study®, there is only one 
Frank constant when the two-fold symmetry of the nematic 
is enlarged to the six-fold symmetry of a hexatic. 

Equilibrium configurations of the director field are found 
by minimizing ^ subject to |n| = 1. This procedure leads to 
h = 0 , where h is the transverse part of the molecular field 
H = —5Fq[/ 5n; h = H — nn • H. Near equilibrium, th e fluid 
stress corresponding to the elastic free energy ^ is then®® 

A 1 

(^ik = -^kldiHi - - {riihk + Tikhi) + - {riihk - rikht ), ( 3 ) 

where Ily^/ = dFQi/d{dkni). In equilibrium, the condition for 
the balance of director torques h = 0 implies the balance of 
elastic forces, —dtp^q -f djojj = 0 , provided the pressure is 
given by peq = —The “tumbling parameter” A is not 


a dissipative coefficient, but is related to the degree of order 
and the type of nematic, with calamitic phases (composed of 
rod-like molecules) tending to have A > 0 , and discotic phases 
(composed of disk-like molecules) tending to have A < 0. The 
value of this parameter further classifies nematic fiuids as ei¬ 
ther “tumbling” (A < 1 ) or “shear-aligning” (A > 1 ). In a sim¬ 
ple shear flow, tumbling nematics continuously rotate whereas 
shear-aligning nematics tend to align themselves at a certain 
fixed angle relative to the principal direction of shear. In 
DSCG, a lyotropic chromonic liquid crystal commonly used 
in experiments on swimming microorganisms in liquid crys¬ 
tals, the tumbling parameter A is a function of temperature 
and has a range A = 0.6 — 0.9®. For comparison, the hexatic 
phase has A = 0 and therefore lacks any of these distinctions. 

Alternatively, we could have formulated the stresses us¬ 
ing the Ericksen-Leslie approach; the connection between the 
Ericksen-Lesl ie ap proach and that used here is discussed in 
the references®®. 

2.2 Governing equations 

The swimming body is modeled as an infinite sheet undergo¬ 
ing a prescribed transverse sinusoidal undulation of the form 
T* = {s/q) sm{qx — cot), measured in the frame moving with 
the swimmer. Here £ is the dimensionless amplitude for the 
swimmer. We focus on transverse waves in the body of this ar¬ 
ticle, but we briefly treat longitudinal waves in the appendix. 

At zero Reynolds number, conservation of mass of an in¬ 
compressible fiuid results in a divergence-free velocity field, 

V • Y = 0, and conservation of momentum is expressed as force 
balance, 

— dip + dj (^afj =0. (4) 

Torque balance is expressed by®® 

aA^,- + (v-V)^,--i[(Vx v) xn]. 

~ ^ i^ij j) ^jk^k F hi/ 1 ( 5 ) 

where 7 * is a rotational or twist viscosit}0 In DSCG, r*/Ai 
ranges from ^ 5 to ^ 5(P. The viscous torque arising from 
the rotation of the director relative to the local fiuid rotation 
balances with viscous torque arising through E and elastic 
torque through — h. We work in the rest frame of the swimmer. 

The no-slip velocity boundary condition is applied on the 
swimmer surface, and as y ^ 00 the flow has uniform velocity 

V = f/*x where —U* is the swimming speed. Meanwhile, the 
director field has a preferential angle at the boundary due to 
anchoring conditions. We will study the case of tangential an¬ 
choring at the swimmer surface, with anchoring strength W®. 

2 For comparison with the much simpler hexatic phase, Appendix [a| includes 
the governing equations for a hexatic liquid crystall^. 
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Since we expand in powers of the amplitude, we may write 
this condition to second order in the angle field: 

-Kidye^W{e-d^Y^)=0, ( 6 ) 

where y = describes the swimmer shape, and ([^ is 

evaluated at y = F*. It is convenient to define the dimension¬ 
less anchoring strength w = W/ {qK \). 

Henceforth we treat v, y, and t as dimensionless variables 
by measuring length in units of q~^ and time in units of 
(0~^. The dimensionless viscosities are defined by il\= jll/jl, 
At2 = 7= 7*/M- It is also convenient to introduce 

the wave speed c = colq, which is one in the natural units. The 
ratio of Frank constants is denoted by Kr = KijK^, and we 
define U = t/*/c, and Q = Q"l{(0£^ jq^) for the volumetric 
fiux. The undulating shape of the swimmer takes the nondi- 
mensional form 

(X,F) = (0,e)sin(x-0. (7) 

The elastic response of the fiuid to deformation introduces 
a length-scale-dependent relaxation time, T = 

For small-molecule liquid crystals, typical values are fl « 


10“^ Pas and K^, ~ 10“^^ N. On the length scale of bacte¬ 
rial fiagellar undulations for which q^l /im“^ the relaxation 
time is T ^ 1 ms. Comparing the typical viscous stress Q 
with the typical elastic stress 0, we find the Ericksen num- 
beiE3, written Er = TCO. Note that unlike the Reynolds number, 
which is always small for swimming microorganisms, the Er¬ 
icksen number for a swimming microorganism may be small 
or large. The beat frequencie s and wavenumbers of undulating 
cilia and fiagella vary widelyEED^ and for experiments o n bac- 
teria in liquid crystals the Ericksen number can be smafil^^ESI, 
Er ^ 10“^ or largeE3 Er ^ 10^. 

2.3 Leading order fluid flow 

Eollowing Tayloi^il, we pursue a regular perturbation expan¬ 
sion in the wave amplitude £. The stream function if/ is defined 
byY = V X (xf/z ); this form ensures V • y = 0. The stream func¬ 
tion y/ and the angle field 6 are expanded in powers of £ as 
xf/ = exf/W ^ 0{£^) and 6 = £0^^^ ^£^6^^^ ^0{£^). 

Eorce and torque balance from and § at 0(e) are given 
by 


V4y/(1) 


4Mi 
2 + 1^2 


d^dy 


1 

(2 + ;U 2 )Er 


{( 1 + [Kr{\+x) + \-xyd^e^^'^+Kr{\-x)d^e^^">^ 

a,e<" + (««'■>+W;«'‘>) 


0 , ( 8 ) 

0. (9) 


Equations ^ and ([^ are solved by = Re[i^(^)] and = 
Re[0*^^)], where 


3 



7=1 

3 

(10) 

0 (1) 

II 

+ 

(11) 


Insertion of and ( pTj ) into ^ and ([^ results in a cubic 
equation for m = rj, 

0 = Ao+ mAi + m^A 3 , 

Aq = —2(2 + /i2) + 7[“(I+ ^iEr(2 + /i2)], 

Ai = m{^^[y(l H-A)^ “h 2(2-j-/i2)] + 2(4 + 4/ii-f- 2 /^ 2 ) 

+ 7[1 — — 2iEr(2-|-2/ii-f-/i 2 )]}, 

M = Ao + 47A+2i^^[7(-l+A^) - 4-4 ;LIi - 2 ^ 2 ], 

^3 = ^r[7(~ I+ 4 + 2 / 12 ]. (12) 


The velocity field remains finite as y —> 00 if the roots rj are 


taken with negative real part. The relation between the coeffi¬ 
cients Cj and dj follows from and 0: 


, Err[l+A-(1-A)r2] 

^ 2{\-Krr]-\Bry) 


(13) 


and the coefficients Ci are determined by the boundary condi¬ 
tions at first order in amplitude 

dyy\y=0 = 0, (14) 

-dy'>\y=o = (15) 

|y=o = (I 6 ) 

2.4 Second-order problem 

The equations at second order in £ have many terms and are 
unwieldy. However, they are simplified by averaging over the 
spatial period. Since the forcing is a traveling sinusoidal wave 
depending on space and time through the combination x — t, 
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Fig. 2 (Color online) Dimensionless swimming speed C/*/(ce^) vs Er = zco for a swimmer in a nematic liquid crystal with = jj,, ji^ = jU, 
Kr = 1.2, and A = 0.6, with anchoring strengths w = 0 (blue), w = 0.1 (red), w = l (green), and w = 5 (brown). The rotational viscosity 7*/jU 
is 5 (left), 25 (center), and 50 (right). The horizontal axis is rescaled in each plot for enhanced resolution. 


averaging over x causes derivatives with respect to x and with 
respect to t to vanish. Denoting the spatial average by (•) , we 
find 


+ ^ /. (17, 

= g, ( 18 ) 

where / and g are given by 

/ = |( 4 s<'>^ e<'i) + (7e('i. (19) 

s = (v<'> ^ V 9 <'>) - 2 i( 49 ('lv'‘>) - ( 20 ) 


(T) 

for {vx ). The swimming speed U is given by the flow speed 

( [^ at y = oo; 

U = {v^^\=o-aj^[y{l-X)g + {2 + Ai2) J/] dy, (25) 


where to obtain (^) we have integrated by parts. Appendix]^ 
discusses some of the details of calculating this integral. 

We will also be interested in another observable. Unlike in 
the case of an unconfined Taylor swimmer in a Newtonian^^or 
Oldroyd-B fluid at zero Reynolds numbei® there is a net flux 
of fiuid pumped by a swimmer in a liquid crystal. In the lab 
frame, the average fiux is given by 


Q= ( {vx-U)dy^ [ {vP-U)dy-{ysvi^^) 
Jy^ Jo 


y=0 


(26) 


with ki = [Kr{l + A) + 1 — A]/(2 + 112 ) and k 2 = Ky — l. Ex¬ 
panding the no-slip boundary condition to second order, we 
find 

(vi"y,=o = -(l'4vy)|,=o, (21) 

where Y is given by Eq. 0. The second-order part of the 
anchoring condition takes the form 


Note that the second term of Eqn. (26) vanishes for a trans 


verse wave since | 3 ;=o = 0- 


(The second term also van¬ 
ishes for a longitudinal wave, since ys = 0—see Appendix C.) 
Therefore, the fiux is also given to second-order accuracy by 


1 ( 2 ) 


■- ({J?)-U)dy. ill) 


where 


-(aye( 2 ))+w(e( 2 )) 




( 22 ) 


(23) 


The swimming speed and velocity field at second order are 
given by solving (Tf) and subject to the no-slip boundary 

condition and no flow at infinity. The result is 


= {JJ^)\y=o-a ['Y{l-^)g + {2 + ^i 2 )F]dy', (24) 


where F{y) = f{y')dy' and a = 2[7(1 - A)2+2(2+;U2)] K 
The boundary conditions on do not enter the expression 


Note our sign convention: a positive U corresponds to swim¬ 
ming towards the left, opposite the direction of wave propa¬ 
gation (see Eig. [H), while a positive Q corresponds to fiuid 
swept to the right, along the direction of wave propagation. 

3 Results 

3.1 Dependence on Ericksen number 

Eor general Ericksen number the solutions of the governing 
equations to second order do not result in elegant expressions, 
but the swimming speed and fiux are readily found and plot¬ 
ted. The method of solution is described in Appendix In 
the following, we use material parameters that closely mir¬ 
ror the properties of disodium cromolyn glycate (DSCG), in 
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Fig. 3 (Color online) Dimensionless swimming speed U*/(ce^) vs. y*/jU for a swimmer in a nematic fluid with jilf = /t, = n, Kr= 1.2, 

and A = 0.75. The three panels correspond to Er = 0.01 (left), Er = 1.00 (middle), and Er = 100. (right). The colors denote anchoring 
strengths: w = 0 (blue), w = 0.1 (red), w = 1.0 (green), and w = 5.0 (brown). Note that the scale for velocity in the left panel is expanded 
relative to the scales in the middle and right panels. 


which ex periments o n swimmers in liquid crystals have been 
performed^^l^^inEo] choose fli = {12 = I, Kr = 1.2, and 
study the swimming speed as a function of Ericksen number 
Er, anchoring strength w, rotational viscosity 7 , and tumbling 
parameter A. 

Eigure shows the swimming speed as a function of Er = 
TCO [recall T = /i/(i^ 3 ^^)] and 7 *//i = 5, 25, and 50. The 
range of 7 *//i is between ^ 5 and ^ 50 for DSCG®1. Observe 
that when the anchoring strength is weak, the swimming speed 
decreases with Er. This behavior was also seen in the case of 
a swimmer in a hexatic liquid crystal®!. When the anchoring 
strength is strong, w ^ 5, the swimming speed is weakly de¬ 
pendent on Ericksen number, becoming independent of Erick¬ 
sen number for large rotational viscosity (Eig.|^ right panel). 
This weak dependence on Er is suggested by the fact that the 
rotational viscosity enters the governing equations in the com¬ 
bination ( 7 Er)“^ (Eqs.[9 18 20). However, when the rotation 
viscosity is in the low range for DSCG, 7 *//i = 5, the swim¬ 
ming speed increases with Ericksen number when the anchor¬ 
ing strength is moderately strong, w = 5 (Eig. left panel): 
when anchoring is important, the swimming speed increases 
when viscous effects dominate as long as the rotational vis¬ 
cosity is sufficiently low. The increase in swimming speed 
with Er at strong anchoring and modest 7 * is not seen in the 
hexatic liquid crystal®. 


All three panels of Eig. [^indicate that the swimming speed 
becomes independent of anchoring strength when the Erick¬ 
sen number is sufficiently large. Once again, because the ro¬ 
tational viscosity enters always in the combination ( 7 Er)“^ 
the value of Er for which the anchoring strength becomes ir¬ 
relevant is inversely proportional to the rotational viscosity. 
In contrast with the case of a transverse-wave swimmer in a 
hexatic liquid crystal, the swimming speed has a weak but no¬ 
ticeable dependence on Ericksen number when the anchoring 
strength is strong and the rotational viscosity is not too large 
(Fig.@ left panel). But because this dependence is weak we 


can say that the large Er limit is the same as the strong anchor¬ 
ing strength limit. Note that the large Ericksen number limit 
is reached at relatively small values of the Ericksen number; 
in all the panels of Eig. the large Er asymptotic value is 
reached or nearly reached when Er = 1. 

As in the case of a hexatic liquid crystal®, the large Er limit 
is singular, since terms with the highest derivatives in the gov¬ 
erning equations vanish in this limit [See Eqs. ([ 8 }|^, (17 -T^]. 
When elastic stresses are small compared with the viscous 
stresses, it is natural to set the Ericksen number to infinity, 
or equivalently, drop all terms involving the Erank elastic con¬ 
stants. The resulting limiting model is known as Ericksen’s 
transversely isotropic fiuid®. However, this limit is singu¬ 
lar, and therefore Ericksen’s transversely isotropic fiuid does 
not give physical results for the swimming speed. In partic¬ 
ular, Ericksen’s transversely isotropic fiuid would incorrectly 
predict that the swimming speed is independent of rotational 
viscosity at large Ericksen number. Examining the right end 
of each of the panels in Eig. shows that speed depends on 7 
at large Er; we now turn our attention to this dependence. 


3.2 Dependence on rotational viscosity 7 * 

Eigure shows the swimming speed as a function of dimen¬ 
sionless rotational viscosity 7 = 7 * //f for various Ericksen 
numbers and anchoring strengths for a nematic liquid crys¬ 
tal. Eirst note that in accord with discussion of Eig. the 
swimming speed becomes independent of anchoring strength 
for large Er. The middle panel of Eig. shows that the de¬ 
pendence on anchoring strength is very weak even for Er = 1 
(which of course can also be observed from Eig. [^. Sec¬ 
ondly, the swimmer can reverse direction. When 7 * //i is large 
enough, t/* becomes negative, meaning the swimmer swims 
in the same direction as its propagating wave. These qualita¬ 
tive features were also observed in the model for a swimmer 
in a hexatic liquid crystal®. However, Eig. also reveals im¬ 
portant differences between swimming in a hexatic and a ne- 
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Fig. 4 (Color online) Swimming speed in /ims“^ vs 7 *//i for a 
swimmer in a nematic liquid crystal with fii =112 = 1 . = 1 - 2 , 

A = 0.75, q=l rad/^um, T = 1 ms, and w = 0. The colored curves 
correspond to different beat frequencies: ft) = lOrad/s (blue), 
ft) = 50rad/s (red), ft) = 100 rad/s (green), ft) = 200 rad/s (brown). 
The corresponding values of Er = Tft) are all <C 1. 

matic liquid crystal. First, the hexatic swimming speed is al¬ 
ways bounded by the isotropic Newtonian swimming speed® 
|f/*| < ce^/ 2 , whereas the nematic swimming speed can be 
greater than the Newtonian speed. Second, there is a maxi¬ 
mum in the swimming speed as a function of rotational viscos¬ 
ity, as long as the anchoring strength is low enough. The max¬ 
imum is most apparent at small Er, and is in the region of mea¬ 
sured rotational viscosities for DSCG (Fig|^ left panel). The 
maximum is less apparent at higher Ericksen numbers since 
in that regime, the 7 ^ 0 limit of the speed is only slightly 
smaller than the value of the maximum speed. 

Note also that the swimming speed depends on the anchor¬ 
ing strength in the limit of low rotational viscosity. When 
7 ^ 0 , there is a decoupling between the flow fleld and the di¬ 
rector held because the molecular fleld h vanishes in this limit. 
In the problem of swimming in a hexatic liquid crystal®, 
this decoupling is complete, and the swimming speed is the 
isotropic Newtonian swimming speed® t/* = (0£^/{2q) when 
7^0. However, the decoupling is only partial in the case of a 
nematic liquid crystal, since in that case the anisotropic terms 
in the viscous stress depend on the director conflguration even 
when h = 0. When h = 0, the director fleld at each instant is in 
equilibrium, and since this equilibrium conflguration depends 
on the anchoring strength, the stress and ultimately the swim¬ 
ming speed depend on the anchoring strength. In particular, 
the swimming speed does not go to the isotropic swimming 
speed when 7 = 0 in a nematic liquid crystal. 

It is interesting to plot the swimming speed in physical units 
to make the dependence on beat frequency ft) more apparent 
(Fig.g. When the anchoring strength vanishes and 7 * //i is in 
the experimental range of 10 and 100 for DSCG®, the swim¬ 
ming speed depends only weakly on the beat frequency ft): ah 



A 

Fig. 5 (Color online) Dimensionless swimming speed vs tumbling 
parameter A for various Ericksen numbers, = jj., = 1^, 

7 * = 20Id, and zero anchoring strength. 

four curves in Fig.j^cross in this region. 

3.3 Dependence on tumbling parameter A 

Figure shows how the swimming speed depends on the tum¬ 
bling parameter A for various Ericksen numbers. For ah Er¬ 
icksen numbers we And a peak near A = 1 , which marks the 
transition between tumbling and shear-aligning nematic liq¬ 
uid crystals®. The maximum is at A = 1 for moderate to high 
Er, and shifts to slightly higher A when the Ericksen number 
becomes small. 

As mentioned earlier, the general expressions for speed and 
flux are too complicated to display. However, there is a rela¬ 
tively compact expression of the swimming speed in the limit 
of large Ericksen number, which we And by expanding the 
flrst-order solutions of in a Taylor series in Er“\ 

and then inserting these values into (rf\-^25\ to And 

22[2 + Ati(l+A) + Ai2]-r(l+A)(A-l)2 1 \ 

8+4M2 + 2r(A-l)2 VEr/ 

(28) 

This expression confirms our general observation that the 
large-Ericksen number behavior is independent of the anchor¬ 
ing strength w. It also shows that when the Ericksen number 
is large, the swimming speed becomes independent of the ro¬ 
tational viscosity when A = 1. The speed and flux as a func¬ 
tion of tumbling parameter for various rotational viscosities 
are plotted in Eig.[^ Note again that although ( [^ and Eig.[^ 
are appropriate for large Er, they are applicable even to the 
modest Ericksen numbers describing experimental systems, of 
size 1-10, since the swimming speed reaches its high Er limit 
at a low value of Er. We do not have an explanation for why the 
swimming speed becomes independent of rotational viscosity 
when A = 1, but we offer the following remarks. Eirst, as 
mentioned previously, the transition from tumbling to shear- 
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Fig. 6 (Color online) Dimensionless swimming speed (top) and flux 
(bottom) vs A at Er = 1000 for a swimmer in a nematic liquid 
crystal with /ij" = jU, = /i and rotational viscosities given by 
7 */jU = 5 (blue), 7*/jU = 25 (red), and 7*/jU = 50 (green). 


aligning nematic liquid crystals occurs when A = I® Sec¬ 
ond, the governing equations lose some of the highest deriva¬ 
tive terms when A = 1, indicating singular behavior and the 
existence of boundary layers near the swimmer that are thin 
in the y direction. And finally, examination of the first or¬ 
der solutions for the angle field in this limit reveal that the 
angle field simultaneously obeys the strong anchoring and no¬ 
anchoring boundary conditions; in other words, the directors 
align exactly tangential to the swimmer surface, but experi¬ 
ence no torque. 

We close this section by describing the dependence of speed 
on tumbling parameter and rotational viscosity for small Er- 
icksen number and weak anchoring strength, which is also an 
experimentally relevant regime. For the hexatic liquid crys¬ 
tal*^, it has been calculated that the first-order velocity field 
is identical to that generated by a swimmer in a Newto¬ 
nian fiuid, so that for Er C 1 and w = 0 the speed is identical 
to the speed in a Newtonian fiuid for any rotational viscos¬ 
ity. In an anisotropic fiuid, however, the flow field can differ 
markedly from the Newtonian counterpart even at first order 
in e, which implies that the speed can differ dramatically from 
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Fig. 7 (Color online) Dimensionless swimming speed vs A for a 
swimmer in a nematic liquid crystal with fii = fi2 = Kr = 1.2, 
w = 0, and Er=T(0 = 0.01, with rotational viscosities y* /fi = 5 
(blue), 'f l\i — 2‘S (red), and y* j/J. = 50 (green). 

the Taylor speed^^t/ = ce^/2, as shown in Fig.[^ In the limit 
A —0, the swimming speed is the same as for a swimmer 
in an isotropic Newtonian fiuid. Note however that there is a 
small but nonzero fiux when A = 0, indicating that the fiow 
field differs from the isotropic fiow field. 

3.4 Swimming and pumping using back flow 

To highlight the role of the nematic degree of freedom in our 
problem, we study swimming and pumping via a mechanism 
in which all fiow is generated by a prescribed motion of the 
directors at a fiat non-deformable wall. The coupling of the 
motion of the directors to the fiow, and vice-versa, is known 
as backfiow. We suppose that some external mechanism oscil¬ 
lates the directors along the wall with the form of a traveling 
wave with wavenumber q and frequency CO, such that the (di¬ 
mensionless) boundary conditions at the wall are 

v|,=o = 0 (29) 

(30) 

Thus, the director configuration rather than the shape is pre¬ 
scribed. For brevity, we call the swimmer with prescribed 
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Fig. 8 (Color online) Dimensionless swimming speed vs Ericksen 
number Tft) for a non-deformable swimmer with prescribed director 
oscillation and jUi = jU 2 = 1, = 1-2, and with 7 */^ = 5 (blue), 

7 */;U = 20 (red), and 7 */^ = 50 (green). 


director configuration a ‘non-deformable’ swimmer, and the 
swimmer with prescribed shape a ‘deformable’ swimmer. 

The swimming speed as a function of Ericksen number is 
shown in Fig. A qualitative difference with a deformable 
swimmer is that the direction does not reverse when the rota¬ 
tional viscosity is large; in fact, increasing the rotational vis¬ 
cosity makes the swimmer go faster, as long as the Ericksen 
number is not too large. Given that the large Ericksen number 
limit is singular, we expect a boundary layer in the velocity 
field and angle field when Er is large. In the problem of a de¬ 
formable swimmer, we found that the swimming speed in that 
limit is governed by the strong anchoring condition. Since the 
strong anchoring condition in this case would correspond to 
no motion of the directors at the swimmer surface, we expect 
the speed to vanish as Er increases, as our calculations show 
(Fig-ig. Note also that when 7=0, there is a complete decou¬ 
pling between the fiow field and director field problems, and 
the swimming speed vanishes. 

Figurej^shows the A-dependence of locomotion and pump¬ 
ing for the non-deformable swimmer. The swimmer can swim 
faster than the Taylor swimmer when A ~ 1 and the rotational 
viscosity is sufficiently large. Note that the behavior of the 
swimming speed is qualitatively similar to that induced by a 
swimmer with a deformable body (Fig.|^. The fiux induced 
by the motion of the directors in the non-deformable swim¬ 
mer (Fig|^ is comparable to the fiux induced by a deformable 
swimmer, indicating that at low Ericksen number much of the 
fiux is driven by the backfiow effect. 

4 Discussion and Conclusion 

Because the nematic phase is more anisotropic than the hex- 
atic phase, more parameters are required in its constitutive re- 
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Fig. 9 (Color online) Dimensionless swimming speed (top) and flux 
(bottom) in dimensionless units vs A for non-deformable swimmer 
with prescribed director oscillation and j^i = fi 2 = \, Kr = 1.2, 

Er = 0.2 and with 7 */jU = 5 (blue), 7 * /^u = 20 (red), and 7 */jU = 50 
(green). 


lation. In particular, there are anisotropic viscosities as well as 
different elastic moduli for bend and splay (and twist) director 
configurations. The tumbling parameter A also leads to fur¬ 
ther distinctions, such as tumbling and shear-aligning, which 
do not exist in the hexatic. Therefore, the hexatic model is 
good quantitative approximation for swimming in a nematic 
when the magnitudes of /ii, /i 2 , A, and ^ 1/^3 — 1 are small. 
Except for ^ 1/^3 — 1, these parameters are usually not small 
for the nematic phase. 

To further highlight the quantitative difference in this 
regime, Fig. [T^ shows the difference in speeds between a 
highly calamitic, near-aligning transition nematic fiuid and a 
hexatic for small values of the Ericksen number and a range 
of rotational viscosities. For these parameters, the swimmer 
in the nematic fiuid can travel at much greater speeds than its 
companion in a hexatic fiuid. 

In this work we extended Taylor’s model of an undulat¬ 
ing sheet locomoting by means of small-amplitude traveling 
waves in a Newtonian fiuid to the case where the ambient fiuid 
is a twist-free nematic liquid crystal. By considering coupled 
equations for the local nematic director and velocity fields and 
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Fig. 10 (Color online) Difference in non-dimensional swimming 
speeds versus 7* /jU and Er between a highly calamitic nematic fluid 
(A = 0.9) with jU*/^ = jU*/^ = I, Kr = 1.2, and a hexatic fluid. 
Both fluids are subject to weak anchoring conditions (w = 0). The 
horizontal and vertical scales are linear, not logarithmic. 


expanding perturbatively in the amplitude we were able to de¬ 
rive general formulas for swimming speed and volumetric flux 
induced by the Taylor sheet. 

Many of the surprising qualitative features, such as rever¬ 
sal of swimming direction for high rotational viscosity, the 
presence of non-zero volumetrix flux, and a convergence to a 
strongly-anchored solution for all anchoring strengths at high 
Ericksen number, have also been seen in the case of a hex¬ 
atic liquid crystal filmISl. However, the effects of anisotropy 
tend for general material parameters to enhance the swimming 
sp eed, a s can occur in for swimming in porous or elastic flu- 
idsEmni^ shear thinning fluids®!, or near rigid walls®!. This 
speed augmentation by anisotropy can be pronounced, partic¬ 
ularly in the low-Ericksen number regime. Our results show 
that the distinctive properties nematic liquid crystals, such as 
backflow, can be exploited to develop novel methods of swim¬ 
ming and pumping in anisotropic fluids. 
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A Hexatic equations 

Eor comparison purposes, we include here the governing equa¬ 
tions for the hexatic phase: 


■ Vp+- KV ■ (veve) +1 v x (zv^e) = o, 

1 K 

(9re+v-ve--z-vxv = -v^e. 

2 7 


(31) 

(32) 


B Details of calculating the swimming speed 

The calculation of the swimming speed, which enters at 
O(e^), depends on a cumbersome but straight-forward com¬ 
bination of the flrst-order flow and director flelds via 
The real part of the first order stream-function in ( p^ may be 
written as 


= + , (33) 

where the overbar denotes the complex conjugate. The direc¬ 
tor angle field is similarly treated. In we require such 

quantities as (OxOyy), which we obtain via 


3 _ 

d^eyy= £ X 

j,k=l 

(dkrie''^y+'^^-*^ + . (34) 

The horizontal mean over one period is then given by 

h = {eAy) = \ t 

^ j,k=l 

(35) 

The final integration in ( [25] ) is now easily performed; for ex¬ 
ample, we have 



djdktj 

(r-j + rA. ’ 


(36) 


and the other contributions are deduced in the same fashion. 
The end result is a cumbersome algebraic expression but one 
that is easily evaluated for all parameter values. 
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Fig. 11 (Color online) Dimensionless swimming speed vs. Er for 
various rotational viscosities for a swimmer with a longitudinal 
wave in a nematic liquid crystal with = fi , Kr = 1.2, 

A = 0.6, and zero anchoring strength iv. 


takes a simple form, with the speed exactly as the transverse 
case but with opposite direction: 

TT— + +M 2 ] - 7(1 “ 1)^ I 1 A 

8 + 4ai2 + 27(A-1)2 \ EtJ 

(43) 

The entrained flux in this limit is likewise of same magnitude 
but opposite direction 

The longitudinal swimmer in a nematic is slower than the 
transverse swimmer. However, a longitudinal swimmer in a 
nematic is very different from a longitudinal swimmer in a 
hexatic. In the case of swimming in a hexatic liquid crystal, 
the longitudinal swimmer’s speed differs from the isotropic 
speed by only a few percent®. Figure[^ shows that the differ¬ 
ence between the isotropic and nematic swimming speed can 
be significant, especially at higher values of rotational viscos¬ 
ity. 


C Swimmer with a longitudinal wave 


References 


For completeness, we also present some results for a swimmer 
with a longitudinal waveform, 

(X,T) = (e,0)sin(x-0, (37) 


in dimensionless form. Many of the equations needed for cal¬ 
culating the speed and fiux are the same as in the case of the 
transverse swimmer. Here we list the equations that must be 
modified. The boundary conditions at first order in amplitude 
£ for the longitudinal swimmer are 


T 

00 

1 

II 

0 

II 

(38) 

0 ^ 

II 

0 

II 

1 

(39) 

-^w6^^'>\y=0 = 0. 

(40) 


The boundary condition for the fiow field at second order is 


(vF))i,=o = -(x4vyoi,=o, 


d)\ 


(41) 


with X given by ( [Tt] ). The anchoring condition remains the 
same as Eq. ([22|), with the replacement of S [Eq. ([23))] with 


s = {xdjdye^^'> -wxdxe^^'>) 


y=0 


(42) 


The swimming speed vs Ericksen number is shown in 


Eig. [TT] Eor most values of Er and rotational viscosity, the 
swimming speed is negative, meaning the swimmer moves in 
the same direction as the propagating longitudinal wave, just 
as in the case of a longitudinal swimmer in an isotropic New¬ 
tonian liquid, where the swimming speed is f/* = —ce^/2. As 
in the transverse case, the swimming direction can reverse if 
the rotational viscosity is sufficiently high. There is no simple 
formula for the swimming speed for generic values of the pa¬ 
rameters, but the swimming velocity at high Ericksen number 
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